Time evolution of entangled biatomic states in a cavity 
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We study the time evolution of entangled states of a pair of identical atoms, considered in the 
harmonic approximation, coupled to an environment represented by an infinite set of free oscillators, 
with the whole system confined within a spherical cavity of radius R. Taking the center-of-mass 
and the relative-position coordinates, and using the dressed-state approach, we present the time 
evolution of some quantities measuring the entanglement, for both limits of a very large and a small 
cavity; the chosen examples are simple and illustrate these very distinct behaviors. 
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Introduction 

Interacting or noninteracting parts of a quantum sys- 
tem can share entangled states that hold quantum cor- 
relations [l|, 0. An recent review covering all relevant 
aspects of the subject is in Q . Entanglement is a purely 
quantum phenomenon due to the attribution of physical 
meaning to superposed states, a concept with no cor- 
respondence in classical physics. Entanglement means 
that individual parts of a quantum system are not inde- 
pendent of each other, even if they do not interact, and 
their quantum properties are described by their common 
wavefunction. In particular, entanglement properties of 
bipartite systems have been largely investigated along re- 
cent years. Within this context, we consider here a sim- 
ple biatomic system, in which each atom is modeled by a 
harmonic oscillator. In this case, studies have been per- 
formed for a noninteracting bipartite systern, with differ- 
ent approaches, as for instance in Refs. 0, [Hj and when 
an interaction between the oscillators is taken into ac- 
count 

In this Brief Report, we study the time evolution of a 
superposition of two biatomic states of identical atoms 
which interact indirectly via the coupling with the har- 
monic modes of a field force representing the environ- 
ment. We consider the two atoms in the harmonic ap- 
proximation and assume that the whole system resides 
inside a spherical cavity of radius R. Our basic objects 
will be dressed states, corresponding to the atoms dressed 
by the field. The biatomic system will be consistently de- 
scribed by the pair consisting of the "center-of-mass" and 
the "relative-position" oscillators, a procedure alrea^ 
employed in the literature, as for instance in Refs. Q ■ 
In our case, these oscillators will be appropriately dressed 
by the field. We will consider the entangled state formed 
by the superposition of two kinds of states: one state 
in which the "center-of-mass" oscillator is in its first ex- 
cited level and the "relative-position" oscillator is in the 
ground state; this state is superposed with another state 
in which the oscillators have their role reversed. We will 
be concerned by the system at zero temperature, i.e. all 
the field modes are in their ground states. Actually, it 



has been shown in [l^ [l^ that thermal effects on dressed 
oscillators are important only for high temperatures, they 
are negligible for room temperatures. This means that we 
can make the approximation of taking the time evolution 
of dressed states, as established for zero-temperature Q. 
Thus, our results apply in fact for temperatures up to 
the order of room temperatures. 
The model 

We start from a slightly modified version of the Hamil- 
tonian used in Refs. [3, describing two atoms A and 
B, in the harmonic approximation, coupled to an envi- 
ronment, the whole system being contained in a perfectly 
reflecting cavity of radius R ^7|] , 



H = 



2 I 2 2 

Pa + ^aQa 



N 



1 2 
^B^lB 



N 

E 

fe=i 



[pk 



2 2\ 

^kik) 



\Fi(r\AqA + ■>iBqB)uJkqk, 



(1) 



fe=i 



where the limit N ^ oo will be understood. Since 
the atoms are identical, we write uja = ojb = wo ^-nd 
Va = Vb ^ V ~ 2-\/(7Aw/7r, where g is a constant with 
dimension of frequency measuring the strength of the 
coupling and Aw = nc/R is the interval between two 
neighboring frequencies of the field [l^ . 

For identical atoms, we define new coordinates, g+ 
(center of mass) and q- (relative position) , such that 
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and corresponding formulas for momenta. Then in terms 
of p± , q± , the Hamiltonian is written as H = + H 



where iJ_ = ^ [pt + u^nQ-l and 
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We see that the "center-of-mass" oscillator, q^, couples 
to the field while the "relative-position" one, q^, oscil- 
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lates freely. The Hamiltonian describes a single os- 
cillator linearly coupled to the field. Therefore, for the 
system ((7+ © field), composed by the "center-of-mass" 
oscillator coupled to the field, we can define the renormal- 
ized frequency ui, dressed coordinates and dressed states. 
Dressed coordinates and states are defined by transfor- 
mations of the normal- mode coordinates and of the eigen- 
states of the diagonalized Hamiltonian, which depend on 
the normal-mode frequencies and on the renormalized os- 
cillator frequency. An important aspect is that excited 
dressed states decay with time while normal modes are 
stationary. Details of this formalism are presented in 
Refs. 0, [H, [H and references therein. 
Time evolution of entangled biatomic states 
We shall describe the pair of atoms by the pair center- 
of-mass and relative-position oscillators. Consider the 
product states |r^o"') = |r^+'> ® \t\^^) and |r[,t"^) = 

® Ir^^-*) in which, the dressed center-of-mass os- 
cillator is in the first excited level, while the relative- 
position oscillator is in its ground state, and vice-versa, 
respectively. The dressed first-excited state of the center- 
of-mass oscillator evolves in time according to [1, [l^ 



and we obtain 



rr'(i))-E/+'^W|rr(o)>: 



(4) 



with — 1, where the label v runs over the 

center-of-mass oscillator (+) and the field modes {{k}). 
The quantity f+v{t) thus represents the probability am- 
plitude that the excitation is at the v-fh dressed oscillator 
at time t. On the other hand, the first excited state of 
the relative-position oscillator is stationary. This implies 



that the time evolution of the state |r 
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is governed 

by Eq. (jH), while the state \T)^ ') remains stationary. 
We now consider at t = the family of states 

l*^'^(o)> = ir^-^o)) + v^e'*|r^r^(o)) (5) 

which belongs to the Hilbert space H^^^k] ® H-, rep- 
resenting states of the system of the two atoms coupled 
to the environmental field. At the instant t, the density 
matrix corresponding to Eq ([5|) is written as 

git) = \^^^it)){^^^it)\=c\r[t-\t)){r[t-\t)\ 

+ Vai^) e''^\T^^r^){r[V{t)\ 

+ y!(W)e-*|r(+-\t))(r(r)|. (6) 

We adopt a more explicit notation and write |rj^ (t)) = 



|l+(i),0_;0,0,...) and |r^t 



(+-)\ _ 



|0+,1_;0,0,...). 



To analyze how the two-atom state evolves in time, 
we consider the reduced density matrix, obtained 
by taking the trace over the field modes, p{t) = 
J2kAki,k2,...\git)\h,k2,...). Using Eq. Q, 

\l+{t), 0_; 0, 0, . . .) = /++(t)|l+(0), 0_; 0, 0, . . .) 

+ E/+,(t)|0+,0_;0,...,0,l„0,...), 



where 
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As it should, Trp(i) = a(<) + 6(t) + cit) = 1. The degree 
of impurity of the state ([7]) is given by £> = 1 — Trp^ it) ; 
then, it follows that 

Dii, t) = 2^(1 - - e + m+it)\'). (13) 

To find the time dependence of the above quantities, 
we need to evaluate the function which governs 

the behavior of the system. There are two significantly 
different situations, depending on the size of the cavity. 
We shall analyze the limit of a very large cavity (i? — > 00) 
and the case of a small cavity. In any case, independently 
of emission frequency, we must have < |/++(i)P < 1- 
Large cavity: for a very large cavity, we have [20| 
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+ iGit-uj,g), (14) 



g'^ and the function G is given by 
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We will consider > 0, which includes the weak cou- 
pling regime g^ <^ uP' . 

Small cavity: for a finite cavity, the spectrum of eigen- 
frequencies is discrete, and the continuum language used 
in the case of large cavity is not valid. However, for a 
very small cavity, with a radius R much smaller than the 
coherence length ■kc/ g (i? <C ttc/ g), we can obtain /4.+ (t) 
to a good approximation as [20| 



f++it) 
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where 6 — gR/irc 1 is a dimensionless parameter, 
characterizing the smallness of the cavity. These two 
extreme cases are illustrated in Fig. [TJ where we plot the 
degree of impurity as a function of the time. 
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FIG. 1: Degree of impurity as function of time, Eq. (|13|l . 
taking cj — 1.5 and g = 1.0 (in arbitrary units), for states 
with 5 — 0.5: for a very large cavity (dashed line) and for a 
small cavity with S = 0.1 (full line). 



Entanglement 

Let us now examine the degree of entanglement of the 
states described by the density matrix (O and how it 
evolves in time. Let us initially calculate the concur- 
rence [2l| associated with the density matrix p. The 
"spin-flipped" state is p = ((72 ^ cr2)p*{<^2 ^ (^2), where 
— i \ . 

„ is the Pauli matrix. Taking the basis 



02 = 

{|0+,0_), |0+, |l+,0_), we find 
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where we have used that h{t)c{t) — \d{t)\'^ . The concur- 
rence of a mixed bipartite state p is given by [2l| 



Cp = max{0, Ai — A2 — A3 — A4} 



(18) 



where the A^s are the square roots of the eigenvalues of 
the non-Hermitian matrix pp, written in decreasing or- 
der. In our case, there is only one nonvanishing eigen- 
value of pp, given by \{t) = Ab{t)c{t). We then find the 
concurrence as 



C,(i) = 2Vfo(iWi)=2Ve(W)l/++(i)l- (19) 
Another measure of entanglement is the negativity [2^ . 



|23|, which can be defined by 



IP^ili 
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where p^^ is the partial transpose of the bipartite mixed 
state p and || • ||i denotes the trace norm. The trace 
norm of an operator O is defined by ||0||i = TrVOOt 
which reduces, for hermitian operators, to the sum of the 
absolute values of its eigenvalues. It can be easily shown 
that, for hermitian operators. 
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FIG. 2: Concurrence (full line) and negativity (dashed line) 
as functions of time, for states with ^ = 0.5, u) — 1.5 and 
g = 1.0 (in arbitrary units), in a very large cavity. 
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FIG. 3: Same as in Fig. [H but for a small cavity with 5 = 0.1. 



so that the negativity is given by the absolute value of 
the sum of the negative eigenvalues of p^~ , that is, 



Mp^2 



E^. 



(21) 



Xj<0 

The partial transpose of p, p^^ , is given by 

a{t) d{t) 

b{t) 

c(<) 
d*{t) 
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(20) whose eigenvalues are \i{t) — b{t), X2{t) 



(22) 



c{t), 



X sjt) ^ (ajt) + ^ a^jt) + 4|d(t)P)/2 and Xi{t) = {a{t) - 
\/ a^it) + 4|(i(t)p)/2. We thus obtain the negativity as 



Afpit) = 2|A4(<)| - ^a^{t)+A\d{t)\^-a{t) 



+ (23) 

Notice that, for pure bipartite states, the concurrence 
and the negativity, as defined above, are equal to each 
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other [23|. However, for mixed bipartite states one has 
A/'p < Cp d^] ■ These features are illustrated in Fig. [5] for 
a very large cavity and, in Fig. [3l for a small one, where 
we plot the concurrence and the negativity for some of 
the states (O, as a function of time. 
Conclusions 

Our study of how an entangled non-gaussian biatomic 
state contained in a cavity evolves in time, leads to, as 
an overall conclusion, that the behavior of the system is 
very contrasting in the cases of a very large cavity (free 
space) or of a small cavity. This is essentially due to 
the behavior of the quantity which governs the 

time evolution of the system, in each case. This quantity 
(the probability that the center of mass oscillator remain 
excited at the first level at time i), goes monotonically to 
zero as i — > cx) for a very large cavity, while it oscillates 
indefinitely with a strictly positive minimum for a small 
cavity. These rather distinct behaviors are exhibited by 
the degree of impurity as illustrated in Fig. [1] 

With regard to entanglement, we have examined the 
concurrence and the negativity as measures of entangle- 
ment; these are plotted in Figs. [5] and |3]for a very large 
and a small cavity, respectively. In a very large cavity, 
these quantities decrease monotonically to zero for large 
times and do not show any sudden death. In contrast, 
for a small cavity, both present an oscillatory behavior 
as time evolves and never vanish, having finite values for 
any arbitrarily long elapsed time. This can be related to 
former results in other context, obtained using dressed 



states, for decay and stability of excited atoms in very 
large or small cavities, for both zero or finite tempera- 
ture [ll,[ll,[2i,[2^. 

Nevertheless it is worthy mentioning that, in some 
cases, a finite asymptotic entanglement persists for long 
times in free space, that is in very large cavities, as 
pointed out in Refs. [1, US, IH-Iil- However, in most 
cases, this phenomenon exists when the environment is 
a thermal bath and for initial states that are two-mode 
squeezed (gaussian) states, with a high degree of squeez- 
ing. These situations are different from the case of the 
zero-temperature bath and non-gaussian initial states we 
have considered in this note, for which no asymptotic 
entanglement for large cavities exists. 

In our case, for large cavities the entanglement vanishes 
asymptotically. On the other hand, we have a permanent^ 
oscillating entanglement for small cavities, at all times. 
In the case of a bath as we have adopted here, the ex- 
istence of the entanglement for all times in a sufficiently 
small cavity, in contrast to the fast decay for long times 
in the very large cavity (free space) , favors small cavities 
as candidates to engender phenomena that require en- 
tanglement as a resource. The discussion presented here 
may be extended to other systems relevant to quantum 
computation. 
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